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Limits on phonon information extracted from neutron pair-density functions 
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We explore the possibility of extracting information about lattice dynamics in simple crystal 
structures from the neutron pair-density function (PDF) through inverse data analysis. Contrary 
to the claims by Dimitrov, Louca, and Roder [Phys. Rev. B 60, 6204 (1999)], and in agreement 
with recent work by Reichardt and Pintschovius [Phys. Rev. B 63, 174302 (2001)], we find that 
the PDF alone is not sufficient for constructing accurate phonon dispersions in the entire Brillouin 
zone in systems with complex lattice dynamics. However, our numerical simulations show that for 
monatomic fee and bee crystal structures it is, in principle, possible to obtain phonon moments of 
complex metals as well as phonon frequencies of simple metals within a few percent accuracy. 

PACS numbers: 63.20.-e, 61.12.-q, 61.12.Bt 
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I. INTRODUCTION 

Our goal in this work is to determine whether it is 
possible to obtain high-quality information about the 
phonon dispersion and lattice dynamics over the en- 
tire Brillouin zone (BZ) without having to measure the 
phonon dispersion in single crystals. This is an important 
question for materials where no single crystals are avail- 
able, or when performing measurements at high pressure 
and high temperature. This workj-was motivated by the 
claim of Dimitrov and co-workersEJ that it is possible to 
extract accurate phonon dispersions by properly model- 
ing the phonon system and applying an inverse data anal- 
ysis technique to the measured neutron pair-density func- 
tion (PDF). In principle, the PDF contains all the lattice 
dynamical information, although most of it is lost after 
the integration over dynamical and directional degrees 
of freedom is, performed. Very recently, Reichardt and 
PintschoviusH seriously questioned the results of Dim- 
itrov et al.,u and concluded that the PDF is rather insen- 
sitive to the precise shape of the dispersion curves, and, 
therefore, that there is no hope for extracting lattice vi- 
brations from experimental PDF measurements. SijHiilar 
concerns were raised by Mellergard and McGreevyfl 

Experimental PDF data obtained from neutron pow- 
der diffraction require various corrections to the raw data, 
which introduce additional errors in the analysis. Thus, 
in this paper we studied synthetic data sets to avoid any 
ambiguity in the analysis. pWe tested and utihzed sta: 
dard reverse Monte CarlcQ and Levenberg-Marquard 
methods for the inverse analysis of the PDF spectra to de- 
termine how much lattice dynamical information can be 
recovered. To do that, we generated synthetic PDF data 
sets for various monatomic fee and bcc crystal structures 
from published tables of generalized Born-von Karman 
(BvK) force constants. Q We used a BvK model to de- 
scribe the lattice dynamics since it is easy to implement 
and is well documented in the literature, although it is 
not the best-suited phonon model for metals. This gave 
us full control when testing the robustness and accuracy 
of the inverse analysis methods being used. We empha- 



size that we wished to determine whether a unique re- 
lationship exists between the phonon dispersion and the 
corresponding PDF of the system, not a unique relation- 
ship between the force constants of the phonon model 
and its dispersion or PDF spectrum. 

Our study of the inverse (indirect) problem of extract- 
ing phonons from a given PDF spectrum complement^ 
the previous studies by Reichardt and Pintschovius,cl 
who studied the forward (direct) problem of obtaining 
a unique PDF spectrum from very different phonon dis- 
persion curves. 



II. ANALYSIS 

We calculated the PDF spectrum by Fourier trans- 
forming the|-j:jijgifder-averaged, coherent static structure 
factor, S{q)l 
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dqqsin[qT)[S{q)-l], (1) 



with atomic number density qq. For simplicity, we 
neglected multiphonon processes {i.e., two-phonon and 
higher-order processes) inpthe computation of the dif- 
fuse scattering part of S'(g)llJ Thus, we assumed that the 
one-phonon scattering process is the dominant inelastic 
scattering process in the range of usually measured scat- 
tering vectors, q < J2-jj^35 — 40 A . However, as was 
pointed out earlier ,QEifE3 one needs to be more careful 
about multiphonon contributions when comparing with 
actual experiments. 

We fit the synthetic PDF data sets gsynth(f j) with re- 
spect to the force constants of the model by minimizing 
the function 
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where F is the number of force constants being fitted, 
and N ^ F. The number of spatial points is of the 
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order of iV ~ 1000, and ""(^O = £crQ{'i'*)r,, /ri is an error 
estimate for Qsynthifija We use a relative error ~ 0.03 
at the first PDF peak maximum at distance r*, which 
is typical for experiments with good neutron counting 
statistics and proper background corrections. Values as 
low as £a ~ 0.01 are feasible in high-precision diffraction 
experiments. Since enters in Eq. |2| only as an overall 
scaling factor, none of our results depend on the absolute 
value of Sa- The sums run from just below the first peak 
through r^v — Inm. Extending the sum to r^v = 2nm 
does not lead to any significant changes in our results. 
The model PDF g{ri) depends implicitly on the phonon 
dispersion and, thus, on the fitted force constants. 

One starts the loop of the fitting procedure with a 
(small) set of plausible force constants and calculates the 
phonon frequencies /(ks) and eigenvectors e(ks) on a fine 
mesh in the BZ. The computation of the static structure 
factor in the one-phonon approximation is straightfor- 
ward, once the frequencies and eigenvectors are known. 
Finally, one computes the model PDF by convoluting the 
static structure factor with the instrumental resolution 
function of the diffractometerli3 and compares it with the 
synthetic PDF. An update of the force constants follows. 
These steps are repeated until a termination criterion is 
met, either in the implemented reverse Monte Carlo or 
Levenberg-Marquardt method. Since the Monte Carlo 
method is computationally intensive and slow compared 
to a generalized nonlinear least-squares minimizer, we 
tested the reverse Monte Carlo method only for models 
with a small set of fitting parameters. 

Next we checked the quality of the resulting phonon 
dispersions by comparing the second phonon moment 
/2 and the phonon dispersion /(ks) with the synthetic 
phonon data. The second moment of a monatomic crys- 
tal is defined byO 

with wave vector k, phonon branch index s, and A^bz 
k-points in the summation over the Brillouin zone. 

In order to quantify the goodness of the phonon dis- 
persion we introduce the merit function 

2 _ 1 [/synth(kiSi) — /(kjSi)] 

where Np is the number of phonon frequencies. Their 
estimated errors cTp are taken to be two percent of the 
frequency, but at least 0.04THz. These are typical error 
estimates reported for neutron triple-axis spectroscopy 
experiments. 

III. RESULTS 

We studied in detail elemental Ni, Ag, Al, Ce and Pb 
(fee), as well as Fe and Nb (bcc). Since we obtained 



very similar results for Ni and Ag, we will not discuss 
Ag separately. The results for bcc structures show very 
similar behavior to the ones with fee structure. Ni has a 
simple phonon dispersion typical of monatomic fee crys- 
tals that can be described very well by including only 
the first few nearest-neighbor (NN) shells of interatomic 
force constants in a BvK model. On the other hand, Pb 
shows a complex dispersion (even when neglecting Kohn 
anomalies) that requires long-range forces up through 
8NN shells in a BvK model. Similarly, Fe is the pro- 
totype of a monatomic bee crystal with a simple phonon 
dispersion, whereas Nb possesses a very complex disper- 
sion. 

For each material we used two very different sets of 
initial force constants (set #1 and set #2) to start the 
fit procedure, in order to test the robustness of the mini- 
mization methods. In the fee simulations shown in Figs. ^ 
through H, we initialized the first shell (INN) of force con- 
stants only (i.e., 3 parameters), either by using the cor- 
rect elastic constants (set #1), or by the corresponding 
values of the INN force constants of the synthetic model 
(set #2). Note that in cubic systems the three indepen- 
dent elastic constants are uniquely defined by three force 
constants (see Appendix). In the case of the bee struc- 
tures we initialized the combined four force constants of 
the INN and 2NN shells (with the 2NN force constant 
XY2 = 0), either by using the correct elastic constants 
(set #1), or by the corresponding values of the synthetic 
model (set #2), 

Also, for large numbers of force constants the 
Levenberg-Marquardt algorithm does not always con- 
verge to the global minimum. Instead it gets easily 
trapped in local minima, depending on the initial val- 
ues of the force constants. If the initial values are chosen 
poorly, then this failure is almost unavoidable. 

In the following subsections A through E we address 
in detail the quality and difficulties of fitting the PDF 
curves and investigate the quality of the corresponding 
phonon dispersion curves, phonon moments, elastic con- 
stants, Debye- Waller factors, and PDF peakwidths for a 
large set of elemental materials. We show results for sys- 
tems that range progressively from very simple to very 
complex phonon models by gradually increasing the num- 
ber of nearest-neighbor interatomic forces in a Born- von 
Karman model. In figures ^ through |l^ we present a 
comprehensive study of the inverse problem of extracting 
phonons from PDF curves from elemental cubic materi- 
als. 



A. Goodness of Fit 

It is obvious from our results for Ni, Al, Ce and Pb dis- 
played in Figs. |l| through ||, and for Fe and Nb in Figs. ^ 
through |l2|, that the PDF is rather insensitive to the zone 
boundary phonons. This is best seen for the cases of Al, 
Ce, Pb and Nb, which have complex phonon dispersions 
with anomalies close to the Brillouin zone boundaries. In 
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FIG. 1: Fitting the synthetic PDF of Ni, generated with 
a 4NN BvK force model. Top: PDF fits vs. number of 
nearest-neighbor shells of force constants for two different 
sets (#1 and #2). Center: Relative error of /2 in percent 
of /f = 8.03THZ. Bottom: Fif mre of merit of computed 
phonon dispersions shown in Fig. H 



FIG. 3: Fitting the synthetic PDF of Al, generated with 
an 8NN BvK force model. Top: PDF fits vs. number of 
nearest-neighbor shells of force constants for two different 
sets (#1 and #2). Center: Relative error of /2 in percent 
of /f = 8.28THZ. Bottom: Fij mre of merit of computed 
phonon dispersions shown in Fig. W 
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FIG. 2: (color) Phonon dispersions along high symmetry di- 
rections of the Brillouin zone obtained frorth fitting the PDF 
curve of a generalized 4NN BvK force model,u using fit models 
with INN and up to 4NN shells (set #1). 
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FIG. 4: (color) Phonon dispersions obtained from fitting the 
PDF curve of a generalized 8NN BvK force model,El using fit 
models with INN and up to 4NN and 8NN shells (set #1). 
Only frequencies (circles) with error bars are included in the 
computation of Xphon in Fig. |[ 



other words, the PDF is not sensitive enough to register 
the out-of-phase motions of neighboring atoms. 

Also there is no direct correlation between the good- 
ness of the PDF fit Xpdf ^^^'^ ^^'^ merit function of the 
phonon dispersion Xphon- ^'^^ example, in Fig. ^ it can 
be seen for the metal Pb that Xpdi #1 order 

of magnitude smaller than Xpdf 4NN and 

5NN BvK models. However, their corresponding Xphon 
are reversed. This means that a good PDF fit can result 
in a bad description of the phonon dispersion, and vice 



versa. Similar conclusions follow from our study of Nb in 
Fig. This is a very important finding. Hence it is gen- 
erally not possible to provide a useful mapping between 
PDF spectra and phonon dispersion curves to solve the 
inverse problem. We find that Xpdf ^^'^ Xphon 
related overall only for simple metals (Ni, Ag and Fe) 
and semi-complex metals ( Al and Ce) , and that both fits 
improve asymptotically by adding more force constants 
to the phonon models. This demonstrates that our algo- 
rithm successfully solves the inverse problem for simple 
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FIG. 5: Fitting the synthetic PDF of Ce, generated with 
an 8NN BvK force model. Top: PDF fits vs. number of 
nearest-neighbor shells of force constants for two different 
sets (#1 and #2). Center: Relative error of /2 in percent 
of /j^"* = 2.48THZ. Bottom: Figure of merit of computed 
phonon dispersions shown in Fig. M 



FIG. 7: Fitting the synthetic PDF of Pb, generated with an 
8NN BvK force model. Top: PDF fits vs. number of nearest- 
neighbor shells of force constants for two different sets (#1 
and #2). Center: Relative error of /2 in percent of f^^^^^ = 
1.99THz. Bottom: Figure of merit of the computed phonon 
dispersions shown in Fig. ^ 
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FIG. 6: (color) Phonon dispersions obtained from fttting the 
PDF curve of a generalized 8NN BvK force model,u using fit 
models with INN and up to 4NN and 8NN shells (set #1). 



cases if the PDF spectra can be obtained with arbitrar- 
ily high accuracy. At this point it is not clear if this is a 
robust result that wiU survive once actual experimental 
data sets are studied (by including statistical and sys- 
tematic errors into the computation of S{q)). Thus, it is 
generally not possible to quantify a priori the quality of 
the extracted phonon dispersions based on Xp^f. 

All PDF fits to the synthetic data sets are almost indis- 
tinguishable, as follows from the extremely small values 
of 

Xpdf (see figures). This means that even the crudest 
BvK phonon model with only INN interatomic forces de- 
viates on average less than 0.1 — 1% from the synthetic 
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FIG. 8: (color) Phonon dispersions obtainedr^om fitting the 
PDF curve of a generalized 8NN BvK model,LJ using fit mod- 
els with INN, and up to 4NN and 8NN shells of force con- 
stants (set #1). Only frequencies (circles) with error bars are 
included in the computation of Xphon in Fig. |^. 



PDF data, which is in good agreement with the analysis 
of the forward problem in Ref. |^. In any realistic ex- 
periment, where statistical and systematic measurement 
errors may result in a Xpdf '-^^ order one, it will be nearly 
impossible to attribute a tiny reduction of less than 10""'^ 
in Xpdi to * significantly improved fit. Hence, in a realis- 
tic simulation a PDF fit with a simple INN BvK phonon 
model will be almost indistinguishable from one with a 
more complex 4NN or 8NN BvK model. 
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FIG. 9: Fitting the synthetic PDF of Fe, generated with 
a 5NN BvK force model. Top: PDF fits vs. number of 
nearest-neighbor sheUs of force constants for two different 
sets (#1 and #2). Center: Relative error of /2 in percent 
of /f'"*'' = 8.77THZ. Bottom: Fit fure of merit of computed 
phonon dispersions shown in Fig. hd 



FIG. 11: Fitting the synthetic PDF of Nb, generated with 
an 8NN BvK force model. Top: PDF fits vs. number of 
nearest-neighbor shells of force constants for two different 
sets (#1 and #2). Center: Relative error of /2 in percent 
of /f = 5.86THZ. Bottom: Fij mre of merit of computed 
phonon dispersions shown in Fig. hd 
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FIG. 10: (color) Phonon dispersions obtainei^ from fitting the 
PDF of a generalized 5NN BvK force model,u using fit models 
with up to 2NN and 5NN shells (set #1). Only frequencies 
(circles) with error bars are included in the computation of 
Xphon in Fig. |. 



r H NpN rpH 



bcc Nb 




[OOh] [hhl] [.5.5h] [hhO] 



FIG. 12: (color) Phonon dispersions obtainejd from fitting the 
PDF of a generalized 8NN BvK force model,u using fit models 
with up to 2NN, 7NN and 8NN shells (set #1). Only frequen- 
cies (circles) with error bars are included in the computation 
of Xphon in Fig. 0. 



B. Phonon Moments 

Since phonon moments are an integrated quantity of 
the phonon dispersion curves, it is plausible to expect 
them to be less sensitive to the details of the phonon 
models being used to generate them. This is indeed the 



case, as can be seen for the second moment shown in 
figures ||, ||, 0, ^, and Here we focus on the sec- 
ond moment only, because (1) it is more sensitive to high 
frequencies (near the zone boundaries) than the lower 
moments, and (2) it enters the free energy functional in 
the high-temperature limit, and can be obtained inde- 
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pendently in a specific heat measurement. For simple 
and semi-complex dispersion curves (Ni, Ag, Fe, and Al, 
Ce) the relative errors of the computed second phonon 
moments track the overall goodness of the PDF fit. Un- 
fortunately, this is not true for more complex dispersions 
(Pb and Nb), where we could not establish a correlation 
between the goodness of Xpdf relative error of 

the second moment A/j/Za^"*'' = I/2 - /2^"*''l//2-*'""'- 
However, even systems with complex phonon dispersions 
allow the extraction of phonon moments within a few 
percent of accuracy. Several years ago Knapp et alxB 
arrived at similar conclusions while studying the mean- 
square relative displacement of the central atom in fee 
materials with extended x-ray absorption fine-structure 
measurements. 

Finally, our computations for Ce, Pb, and Nb show 
that the second phonon moment is too insensitive to the 
phonons at the zone boundaries of the Brillouin zone to 
be useful for determining zone boundary phonons. 



C. Elastic Constants 

Dimitrov and co-workerJ3 suggested that by adding 
measured elastic constants as constraints to the PDF 
fit (constrained fit) one can improve upon the extracted 
phonon dispersions. This is contrary to our own results. 
We could not observe any significant changes to our ex- 
tracted phonon dispersions by using constrained PDF 
fits; i.e., constraining the PDF fit to give force constants 
that yield the correct elastic constants does not result in 
better phonon dispersions near the zone boundaries. Of 
course it results in slightly more accurate phonon disper- 
sions near the zone center. This is not surprising since the 
elastic constants determine the long-wavelength limit of 
the phonon dispersions at the center of the Brillouin zone 
(r point) and not at the zone boundaries, where the dis- 
crepancies between the synthetic phonon dispersions and 
the PDF-fitted phonon dispersions are largest. Here the 
extracted phonon dispersions could improve by providing 
additional constraints on zone boundary phonons. 

Furthermore, our analysis shows that, in many cases, 
the unconstrained PDF fits already yield elastic con- 
stants that deviate only a few percent (approximately 
1 — 6%) from the values of the elastic constants of the 
synthetie data sets, when a significantly large number of 
force constants are being used. By "sufficiently large", 
we mean at least a 3NN BvK force model for Ni, Ag, 
Ce and Fe and a 5NN BvK model for Al and Pb. This 
does not work for Nb, however, where the extracted elas- 
tic constants are more than 10% off, even for the best 
fit. For a definition of the elastic constants in terms of 
a BvK force model for fee and hee crystal structures see 
the Appendix. 



D. Debye- Waller Factor 

Reichardt and Pintschovius suggested that it might be 
possible to improve the quality of extracted phonon dis- 
persions by adding constraints to the PDF fit, e.g., ther- 
mal parameters independently measured by a Rietveld 
analysis. The thermal parameters, which are given by 
the exponent of the Debye- Waller factor, e~^^, measure 
the mean-square atomic displacement {u^). For cubic 
crystals the Debye- Waller exponent simplifies toQ 

2W{ci) = {{ci-nf)=q^ul) = \q^u^), (5) 

where Uq is the component of the displacement vector u 
in the direction of the scattering vector q. This result 
may be expressed another way, 

w{c) = -J^ y y ^^il^coth ( M^ff) ) 

. _|_| ^WM, (6) 

with the Debye- Waller spectral function given by 
, , N(uj) f huo \ 

where M is the atomic mass, a;(ks) is the angular fre- 
quency of mode s, and e(ks) is its normalized eigenvec- 
tor. The phonon density of states is normalized so that 

/■OO 

/ dijjN{uj) = 1 , (8) 
Jo 

and N{lo) = for w larger than the maximum phonon 
frequency. The thermal parameters extracted from our 
PDF fits are remarkably insensitive to the specific form of 
the phonon model. Even for the simplest phonon models 
used, the relative error of {v?) is typically less than 1% 
(less than 0.1% for Ni, Ag, Al, and Pb) and less than 
5% for Nb. Thus, measured thermal parameters, whose 
absolute values are known only within 5%, cannot give 
improved PDF fits or better phonon dispersions. 

At this point a caveat is needed for using elastic con- 
stants (derived from long wavelength modes) and a Debye 
phonon model to estimate the Debye- Waller factor, or 
when comparing Debye temperatures with Debye- Waller 
factors measured at high temperatures. Figs. |l^ and |lj 
show very clearly for Ni and Ce that at room temperature 
the Debye- Waller spectral function W(ti;) estimated from 
the long wavelength modes, is only a crude approxima- 
tion. It consistently underestimates the contribution of 
the intermediate frequency region (transverse modes near 
the zone boundary) , and describes the high frequency re- 
gion (longitudinal modes near the zone boundary) only 
on average. Here the Debye- Waller exponents 2W com- 
puted from a Debye phonon model are 20% too small. 
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FIG. 13: A comparison of the Debye- Waller spectral fuiitftions 
of Ni at 300K for a realistic 4NN BvK phonon modeia with 
its corresponding Debye phonon model. 
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FIG. 14: A comparison of the Debye- Waller spectral fuactions 
of Ce at 300K for a reahstic 8NN BvK phonon modeB with 
its corresponding Debye phonon model. 



compared to the ones using more realistic lattice dynam- 
ical models. The other materials studied show similar 
discrepancies, except in the case of Nb where both the 
Debye and lattice model calculations of the thermal pa- 
rameters accidentally agree within 1 percent. Indeed, it 
is well known that the Debye- Waller factor is rather in- 
sensitive to the detailed form of the phonon dispersion, 
and at high temperatures it depends only on a single pa^ 
rameter, namely, the inverse-squared phonon moment .t£l 
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FIG. 15: The spectral function of the first PDF peak of Ni at 
300K. Insert: Si{uj) /N{uj) is nearly constant for frequencies 
< (j;max/2. This frequency range is significantly smaller 
for Pb and Nb and marks the deviation from a Debye-model 
phonon spectrum. 
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FIG. 16: The spectral function of the first PDF peak of Ce at 
300K. Insert: Si{uj) /N{u)) is nearly constant for frequencies 

UJ < (Jniax/4. 



proach to compute the PDF. They showed in the har- 
monic approximation that the PDF is approximately a 
series of Gaussian peaks, each centered at distance 
with width tXi, 



E. PDF Peakwidths 

For studying local atomic structure properties in semi- 
conductors Chung and ThorpeM used a real space ap- 
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hLu{k.s) 
2kBT 



s keBZ 

X [l — cos(k • r 



■e(ks) 



Lu(ks) 



(9) 
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Here is a position vector of an atom in the i-ih NN 
shell measured relative to an atom at the origin, and 
fi — Yi/ri. In the limit k ■ 3> 1 the cosine term os- 
cillates so rapidly across the BZ that its average van- 
ishes and the peakwidths of far-out atoms contain the 
same information as the Debye- Waller exponent, namely, 
linii^oo erf = o"^ = 'iW{q)/q'^. This explains why the 
PDF fits reproduce the Debye- Waller factors so well, and 
that adding thermal parameters cannot provide any ex- 
tra constraints. 

On the other hand, the width of the first PDF peak 
(which is the sum of the auto-correlation and cross- 
correlation parts of the displacement-displacement func- 
tion) exhibits the largest deviations from the long-range 
value (7^ (which is purely the auto-correlation part of the 
displacement-displacement function), due to the short- 
range nature of correlations of lattice vibrations. Note 
that in the real space approach all multiphonon processes 
have been included. Expanding cos(k-ri) in Eq. ||, it fol- 
lows that at high temperatures and for low frequencies 
the spectral function of af is, in first-order approxima- 
tion, dominated by the phonon density of states. We 
define the spectral function Si{u}) of the i-th peak as 

a} = J^J^ dcjS.icj), (10) 
S^iu;) = >V(cc;)(|n-e(ks)|2[l-cos(k-rO])^, (11) 



IV. CONCLUSIONS 

In conclusion, our study shows that one cannot obtain 
accurate phonon dispersions from an inverse analysis of 
the pair-density function, unless the lattice dynamics are 
simple and fully described by a few phonon parameters, 
as in the cases of fee Ni and Ag and bcc Fe, for exam- 
ple. A semi-quantitative picture of the phonon disper- 
sion may be obtained in simple and semi-complex metals, 
but not in metals with complex phonon dispersions. In 
principle, in simple metals phonon frequencies can be ex- 
tracted within a few percent accuracy (~ 2 — 8%) in the 
entire Brillouin zone, whereas in semi-complex and com- 
plex metals such accuracy applies only to a small fraction 
of the Brillouin zone centered around the P-point with 
wave vector k ^ w/Aa. We found numerically that the 
pair-density function provides an overall account of the 
lattice dynamics by yielding phonon moments within a 
few percent accuracy. In other words, a rather simple 
phonon model suffices to describe the dynamics embed- 
ded in powder diffraction data. Neutron or x-ray PDF 
studies play an important role in the studies of the local 
structure of crystals, but cannot provide deeper insight 
into the dynamics of lattice vibrations. A more promising 
approach for extracting phonons from powders or poly- 
crystals may be the analysis .of jtime-of- flight spectra from 
inelastic neutron scatteringO'EJ 



where (. . is a normalized average over phonon modes 
and k-points at fixed frequency, aj(ks) — uj, and Si(ui) = 
for frequencies larger than the maximum phonon 
frequency. For sufficiently small frequencies we get 
(|fi • e(ks)|2[l - cos(k • ri)])^ - fc^a^ ^ {a/cfto'^, with 
the average sound speed c and lattice constant a. For 
temperatures T ^ huj/2kB^ the width of the first peak 
is a rough measure of the integrated phonon density of 



states times , as can be seen in the inserts of Figs. |ij 
and |l|. 

It is not too surprising to find that the widths of the 
first peak and that of very distant peaks, or equiva- 
lently the Debye- Waller factor, are crudely approximated 
by long wavelength phonons, although for quite differ- 
ent physical reasons. Since for purely elastic neutron 
scattering |-the peakwidths are identical, af = ct^ = 
AW{q)/q^^ and in most metals inelastic scattering (one- 
phonon, two-phonon, and higher-order processes) leads 
only to small corrections Atri of order 10 — 30% in the 
width of the first peak cti, one does not expect large 
variations in <Ti. Hence the magnitude of the peakwidths 
will be approximated in leading order by the long wave- 
length phonons. For example, in the case of Ni we find 
a difference of Acti/cti k, 13% between the fitted PDF 
curve that includes only elastic processes versus elas- 
tic plus one-phonon processes. For Al and Pb we ob- 
tain A(Ti/(Ti w 25%. These one-phonon corrections are 
largest for thjeJirst PDF peak and much smaller for the 
other peaks.l3't3 
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APPENDIX: ELASTIC CONSTANTS IN A 
BORN- VON KARMAN FORCE MODEL 

In a cubic crjpial there are only three independent 
elastic constants.Ej For fee crystal structures with lattice 
constant a these are related to the force constants (up 
through 8NN interatomic shells) by 

aCii = 4A:a:i + AXX2 + IQXX3, + 8YY^ + 16A:A4+ 

36X^5 + AYY^ + 16a:a:6 + ^2XX^+ 

ZTYY-j + %ZZ^ ^ IGX ATg , (A. 1) 

aC44 = 2XXi + 2ZZx + AYY2 + 4X^3 + 20rr3+ 

+ 8ZZ4 + 2xx^ + isrys + 20ZZ5+ 
lexAg 20XX7 -f 4oyr7+ 

hlZZ-j + leFYg , (A.2) 
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a(Ci2 + C44) = 4XFi + SyZa + 32XZ3 + I6XY4+ 
2AXY5 + 32YZe + 96X^7+ 

+ 32YZ7 , (A.3) 

and for bcc structures these are given by 

aCii = 2XXi + 2XX2 + 8XX3 + I8XX4 + 4YYi+ 
8XX5 + 8XXe + 36XXr + 2ZZ7+ 

32XX8 + 8YYs , (A.4) 

aCii = 2XXi + 2YY2 + iXXs + iZZs + 2XX4+ 
20YYi + 8XX5 + 8YYe + 

18ZZ7 + 4x^8 + lerig + 2ozzs , (a.s) 



a(Ci2 + C44) = 
16X^5 



AXYi + 
24:YZ7 



8XY3 + AYZ4 + 24:XY4+ 
f 36XY7 + 32XYs , (A.6) 



where we fohowed the derivation of Squires.E3'til Our elas- 
tic constants agree with those in Ref. ^ where compar- 
ison is possible, except for the term AYY4 in Cn for bcc 
structures. The generalized BvK force matrix of the n-th 
NN interatomic shell is defined by 



(XXfi XYji XZji 
XYn YYn YZ^ 
XZn YZn ZZn 

Its symmetry properties are listed in Table ||. 
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shell 


fee 


bcc 




INN 


[110] 


[111] 






XXi = YYi 


XXi = YYi = ZZi 






XZi = YZi = 


XYi = XZi = YZi 




2NN 


[200] 


[200] 






YY2 = ZZ2 


YY2 = ZZ2 






XY2 = XZ2 = YZ2 = 


XY2 = XZ2 = YZ2 = 





3NN 


[211] 


[220] 






YYs = ZZ3 


XX3 = YY3 






XY3 = XZ3 


XZ3 = YZs = 




4NN 


[220] 


[311] 






XX4 = YY4 


YYi = ZZi 






XZi = YZ4 = 


XYi = XZi 




5NN 


[310] 


[222] 






= = 


XX5 = FYs = ZZs 








XYs = XZs = YZs 




6NN 


[222] 


[400] 






XXe = YY6 = ZZa 


YYe = ZZe 






XYe = XZe = YZ(. 


XYe = XZe = YZs = 





7NN 


[321] 


[331] 








XX7 = YY7, XZ7 = YZ7 


8NN 


[400] 


[420] 






YYs = ZZs 


XZg = YZa = 






XYs = XZs = YZs = 
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